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Abstract 

■ In this paper we study a model of synchronization process on scale free networks with degree- 

degree correlations. This model was already studied on this kind of networks without correlations 
by Pastore y Piontti et al, Phys. Rev. E 76, 046117 (2007). Here, we study the effects of the 
degree-degree correlation on the behavior of the load fluctuations Ws in the steady state. We 
J>^- found that for assortative networks there exist a specific correlation where the system is optimal 

synchronized. In addition, we found that close to this optimally value the fluctuations does not 
depend on the system size and therefore the system becomes fully scalable. This result could 
be very important for some technological applications. On the other hand, far from the optimal 
correlation, Ws scales logarithmically with the system size. 

PACS numbers: 89.75.-k;89.20.Wt;89.75.Da 
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I. INTRODUCTION 



In the last decades the study of complex networks received much attention because many 
real processes work over these kind of structures. Historically, the research was mainly 



ogy affects processes such as epidemic spreadings traffic flow 

6|. Many real networks have 



bcused on how the topo 

2|, |3| , cascading failures and synchronization problems 
structures characterized by a degree distribution P{k) ~ /c"'^ known as scale free (SF), where 
k is the degree or number of connections that a node can have and kmax ^ k > kmin, where 
kmax is the maximum degree, kmin the minimum degree and A measure the broadness of 
the distribution Ql- In synchronization process it is customary to study the fluctuations 
W = |l/A^ ~ (^))^} ^ of some scalar field h, where hi with i = 1,N represent 
the scalar field on node i, (h) is the mean value, is the system size and {.} denotes an 
average over network configurations. These kind of problems are very important in many real 
situations such as supply-chain networks based on electronic transactions {sj], brain networks 

9] and networks of coupled populations in correlated epidemic outbreaks [l^. Pastore y 



Piontti et. al 



ll\ studied a model of surface relaxation with non- conservative noise that 



allows to balance the load and reduce the fluctuations (synchronize) of the scalar fields on 
SF networks without degree correlation. However real networks are correlated in nature, 
and there should be a reason for this feature. One reason could be to enhance some process 
such as the transport and the synchronization through them. The degree-degree correlation 
of a network can be measured using the Pearson's coefficient given by 121] 



M''TltA{3r + ki) 



2 ' 



(1) 



+ k^) - [M-1 + h 

where M is the number of edges of the network and ji and ki are the degree of the nodes of 
the edge i. This coefficient only can takes values in the interval [—1,1]: if r < the network 
is called disassortative (nodes with low degree tend to connect with highly connected nodes) 
while for r > the network is called assortative (nodes tend to connect with others with 
the similar degrees). When r = the network is uncorrelated. As observed in many other 
works the degree-degree correlation affects considerably the processes that occur on top of 
them (l3|-ll5|. 

In this paper we study the effects of the degree-degree correlation on the behavior of the 
fluctuations in the steady state Wg of SF correlated networks with A < 3 for the model 



of surface relaxation to the minimum (SRM) [16] used by Pastore y Piontti et. al [ll|] in 
uncorrelated networks. To study the fluctuations we map the process with a problem of a 
non-equilibrium surface growth 17|, where the scalar field hi = hi{t) represents the interface 
height at each node i at time t. We found that for every A < 3 there exist a value of the 
correlation for which the fluctuations are minimized, i.e, that optimizes the synchronization. 
Close to and at the "optimal" correlation the fluctuations does not depend on A^, but for 
other correlations the fluctuations diverges logarithmically with A^. 



II. MODEL AND SIMULATION 



To construct the networks we use the configurational model (C M) |18[| with a degree cutoff 
kmax = N^^"^ for A < 3 in order to uncorrelate the original network 19|. Then, we choose two 
links at random and with probability p we connect the nodes with higher degree between 
them and the two with smaller degree to each other to obtain r > 0. For r < 0, we connect 
with probability p the node with highest degree with the one with lowest degree and the 
other two between them. In both cases we do not allow self loops or multiple connections. It 
is known that algorithms that generate clustering (the probability that two connected nodes 
have another neighbor in common) produce degree-degree correlation, but the algorithm 
used here produce degree-degree correlation without introducing clustering 2^. In this 
way, we can study the effects of the degree-degree correlation on SF networks isolating them 
from clustering effects. A side effect of this algorithm is that for SF networks the range 
of Pearson's coefficient that can be generated cannot span the total domain r G [—1,1]. 
Nevertheless, the range that can be obtained is enough to observe how change the scaling 
of the fluctuations with the system size when correlations are introduced. For all the results 
in this work we use kmin = 2 in order to ensure that the network is fully connected 2l| . We 
present the results for A = 2.5 but we checked that for 2 < A < 3 they are qualitatively the 
same. The reason to investigate only 2 < A < 3 is because almost all the real SF networks 



fall in this range of va^ 
In the SRM model 



ues o: 



f A. 
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16| , at each time step a node i is chosen to evolve with probability 
1/A^. Then, if we denote by Vi the nearest neighbor nodes of i, the growing rules are: (1) 
if hi < hj Vj E Vi ^ hi = hi + 1, else (2) if hj < /i„ Vn 7^ j G f j ^ hj = hj + 1. For 
the simulations we start with an initial configuration of {hi} randomly distributed in the 



3 



interval [0, 1]. 

In Fig. [T]we show, in log-linear scale, Ws as a function of for different values of r for 
A = 2.5. We can see that for some values of r, Ws has a logarithmic divergence with 



while for other values of r, Ws does not depend or has a weakly dependence on A^. 
change of behavior means that the scaling of the fluctuations not only depends on A 



'his 



111, 



it also depends on the correlation of the network. These results are in agreement with Ref. 



ll| , where for uncorrelated, or slightly disassortative networks, Ws scales as In for A < 3 
(r = —0.05 in Fig. [1]). Notice that the relation between r and p has finite size effects (See 
Fig. [2]). For this reason if we want to fix r we must select different values of p for each 
system size. 

In Fig. [3] we plot W^ as a function of r for A^ = 5000. Each data point was obtained from 
the linear fitting of W'^{t) in the saturated regime for each r value for 3000 realizations, 
a task of very time consuming. We can see that there is a positive value r = r^m that 
minimizes (optimizes) the fluctuations. In the inset figure we show Tmin as a function of A^. 
We can see that for large system size (A^ ^ 3000) the optimal correlation is independent of 
A^. The dashed line represent the linear fitting of Tmin for large A^, from where we found 
that rmin ~ 0.335 for A = 2.5. This means that for the optimal correlation the fluctuations 
in the steady state do not depend on the system size. This is an important result because 
for r close to r^nin the whole system is scalable with A^. As an example, suppose that we 
have a cluster of computers connected as a SF network with r ^ Tmin, and that the excess 
of load of the cluster of computers is sent to the first neighbors in one time step as in our 
model. In the optimal correlation, as we show, the fluctuations are independent of A^, so we 
could increase the number of computers in our system as much as we want without losing 
its synchronization. 

In order to explain this behavior we compute the local contribution to the fluctuations 
due to all nodes with degree k in the steady state, given by 

1 ^ 
W? = = V W^ 

and therefore the total fluctuation can be computed as 

^max 

W^= E P{k)W^. (2) 
In Fig. Hla) we show W^ as a function of A; for A = 2.5, A^ = 5000 and different values of r. 



We can see that for nodes with high degree, decreases as r increases. This is due to the 
fact that when r increases nodes tend to connect with others with similar connectivities and 
since the high degree nodes are few and tightly packed, the one step relaxation is enough for 
even out all their heights, balancing better the load and enhancing the synchronization. On 
the other hand, for low degree nodes we observe that has a minimum for the optimal 
correlation, as shown in Fig. H] b). In SF networks the majority of the nodes have low 
connectivities and if r > the average distance between those nodes becomes bigger that 
one 13|. As the relaxation is only to first neighbors, different parts of those chains will 
have very different heights. As a consequence, for some values of r > not all the low 
degree nodes will be completely synchronized between them. For this reason the optimal 
correlation is positive and smaller than one since low degree nodes are connected to some 
nodes with high degree allowing to speed up the relaxation and smoothing out the interface 
among them. 

In order to prove this, in Fig. [5] we plot the cumulative 

k 

w^{k)= ^y: w',,p{k'). (3) 

k =kmin 

As we can see from the plot, as r increases the contribution to Wg of high degree nodes 
decreases, being the main contribution to the fluctuations due to low degree nodes for 
positive r and the smaller for r^m. This is why the global fluctuations is minimal in the 
optimal correlation. Also from the same plot we can understand why close to the optimal 
correlation the system does not depend on N . For k > k*, Eq. ([3]) can be rewritten as 

k 

w^{k) = w^{k*)+ E w',,p{k'), 

k'=k*+l 

where k* is the upper value of k that separate two different regimes for W^{k). The second 
term can be replaced by A{k), where 

A ^-1 \ In fc, r < Tjnin and r > r^jn j , 
A{k) ~ (4) 

const., r ~ Tmin ■ 

Then Eq. (jj]) is given by 

= W^,{kma.) , 

where kmax = „ . „ 

In y/N ~ In A^, r < and r > Vmin ] 

(5) 




The logarithmic divergence far from the optimal correlation is a consequence of that the 
contribution of high degree nodes cannot be disregarded. However in r = r^m only the low 
degree nodes contribute to the global fluctuations of the system allowing the scalability with 
the system size. 

III. SUMMARY 

In this paper we study the effects of degree degree correlations on the behavior of the 
fluctuations Wg for the SRM model in SF networks with 2 < A < 3. We found that there 
exist an optimal value of the Pearson's coefficient < Vmin < 1 (assortative networks) 
where the system is optimal synchronized. We also found that for values close to r^m the 
fluctuations does not depend on N, i. e, it is scalable. Moreover for r < Vmin and r > Tmin the 
fluctuations diverge as a logarithmically with the system size N. Then the scaling behavior 
of Ws with the system size depend strongly on the correlation of the network. 

The optimal synchronization found for assortative networks is in our model a topological 
effect due to the correlations. This is an unexpected result because in many researches it 
was found that disassortative networks (communication networks) are better for transport 



and to synchronize oscillators 15 1. 
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FIG. 1: Ws as a function of N in log-linear scale for A = 2.5 and kmin = 2 for r = —0.237 {Q), 
—0.05 (□), 0.103 (o), rmin = 0.335 (A) and 0.386 (*). The dashed lines represent a logarithmic 
fitting and the dotted lines a linear fitting. 
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FIG. 2: r as a function of p for A = 2.5 and kmin = 2 for AT = 1000 (straight line), 3000 (dotted 
line), 5000 (dashed line) and 7000 (dot-dashed line). The straight horizontal line with the arrow 
indicate the values of p used for r = 0.2. 
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FIG. 3: Wg as a function of r for A = 2.5 and N = 5000. The arrow indicate the position rmin- In 
the inset figures we plot Vmin as a function of N in symbols. The dashed line represent the linear 
fitting in the region where rmin is independent of N: rmin = 0.335 for A = 2.5. The averages were 
done over 3000 realizations. 
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FIG. 4: vs for A = 2.5, r = -0.255 (o), -7.10-°^ (□), 0.335 (rmin) {<>) and 0.386 (A) b) An 
amplification of a) for low degree nodes. The dashed lines are used as guides. All these results are 
for TV = 5000 and 5000 realizations of the networks. 
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FIG. 5: Log-linear plot of W^ik) as a function of for A = 2.5, r = -0.255 (o), -T.IQ-^^ (□) and 
0.335 {vmin) (o)- We can see that for k > k* , where k* 10 the asymptotic behavior of W'^{k) 
goes as In A; for r < Vmin and r > Vmin a^nd goes as a const, for r 2± Vmin- 
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